Harrison-Zel'dovich primordial spectrum is consistent with observations 
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Inflation predicts primordial scalar perturbations with a nearly scale-invariant spectrum and a 
spectral index approximately unity (the Harrison-Zel'dovich (HZ) spectrum). The first important 
step for inflationary cosmology is to check the consistency of the HZ primordial spectrum with 
current observations. Recent analyses have claimed that a HZ primordial spectrum is excluded at 
more than 99% c.l. Here we show that the HZ spectrum is only marginally disfavored if one considers 
a more general reionization scenario. Data from the Planck mission will settle the issue. 

PACS numbers: 98.80.-k 95.85. Sz, 98.70.Vc, 98.80.Cq 
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One usually models the dynamics of primordial infla- 
tion by the dynamics of a scalar field, know as the in- 
flaton, evolving under the influence of a scalar potential. 
The scalar and tensor perturbation spectra produced by 
inflation depend upon the value of the inflaton potential 
during inflation and how rapidly the scalar field evolves 
during inflation. Observational information of the spec- 
tra can be translated into information about the inflaton 
potential, thus about physics far beyond the scales of the 
standard model of particle physics 

Cosmic microwave background (CMB) observables 
that map onto the the potential and its slope are 1) n, 
the spectral index of comoving scalar perturbations P-jz: 
n—1 — d\n P-ji/dlnk, and 2) r, the ratio of tensor (grav- 
itational wave) perturbations to Pn [2] . 

CMB temperature and polarization anisotropics are 
the most promising route to explore the physics behind 
inflation. Recent measurements by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) satellite seven-year 
mission [3] , combined with the ground-based and balloon- 
borne experiments such as BOOMERanG [4], QUAD 5], 
ACBAR and BICEP have sharpened our knowl- 
edge of some key inflationary parameters. 

With regard to the dynamics of inflation, a hotly de- 
bated question is whether the case of n = 1 is signiflcantly 
excluded by current observations (see e.g., Refs. (3|-[l3|). 

The scalar spectral index, n, has been recently con- 
strained to the value n = 0.963 ± 0.014 at 68% c.l. by 
the WMAP seven-year dataset (WMAP7) [1], disfavor- 
ing the value of n = 1 at about two standard deviations. 
A combined analysis with galaxy clustering data gives 
n ~ 0.963 ± 0.012 at 68% c.l., rulin g ou t n = 1 at more 
than three standard deviations ([sj. [l3|). 

Compelling evidence that n 1 would be quite reveal- 
ing for two reasons. First, it would provide an indication 
for the dynamical evolution of the expansion rate of the 



universe as perturbations are being produced. Secondly, 
in many models of slow-roll inflation the amplitude of the 
tensor perturbations are proportional to [jt.— 1|. Thus, the 
larger the departure of n from unity, the more likely ten- 
sor perturbations would be within observational reach. 

As it is well known, a scalar spectral index with an 
exact value of n = 1 corresponds to the phenomenologi- 
cal model proposed by Harrison, Zel'dovich, and Peebles 
14|.-^ While one can construct inflationary models that 
give n = 1 either approximately [3] or exactly [l^ , they 
are less than compelling. In fact, observations pointing 
to n exactly unity may indicate that the origin of cosmic 
perturbations lies in some unknown fundamental process 
and may not arise from inflation. 

In this paper we show that the value of n, and its uncer- 
tainty, derived from CMB datasets are sensitive to how 
one models the reionization of the universe. In particu- 
lar, the statement that n = 1 is observationally excluded 
is very much weakened if one treats reionization in a gen- 
eral manner. 

It is well known from a large set of astrophysical ob- 
servables that after primordial recombination (which oc- 
curred at a redshift of z ~ 1100) the universe "reionized" 
at a redshift z > 6. Current constraints from CMB data 
assume a sudden and complete reionization at a redshift 
Zr- This parameter is usually assumed to be in the range 
4 < Zr < 32, and the constraints on n are obtained after 
marginalization over Zr. The electron ionization frac- 
tion Xe{z) is parametrized by Zr such that for z <C z^ 
Xe(z) = 1 {xe{z) = 1.08 for z < 3 in order to take into 
account Helium recombination) and Xe{z) — 2x 10~^ for 



^ For tlie purpose of our paper, we will assume that the tensor 
modes are small enough as to be irrelevant. 
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z > Zr, i.e., joining the value after primordial recom- 
bination with a smooth interpolation. We refer to this 
common procedure in what follows as "sudden" reioniza- 
tion. 

However, the precise details of the reionization process 
are not very well known, and therefore the reionization 
history of the universe at those redshifts could have easily 
been very different. 

As already remarked in Refs. and [3], the assump- 
tion of a more general reionization scheme could affect the 
cosmological constraints on n. It is therefore timely, es- 
pecially in view of the recent CMB and galaxy clustering 
data that have improved the bound on n, to investigate 
the impact of the reionization history on the current con- 
straints on n in a more general reionization scenario. 

We adopt two methods for parametrizing the reioniza- 
tion history. The first method, developed in Ref. [l^l, is 
based on principal components that provide a complete 
basis for describing the effects of reionization on large- 
scale E-m.ode polarization. Following Ref. [13], one can 
parametrize the reionization history as a free function of 
redshift by decomposing Xe{z) into its principal compo- 
nents: 



Xe{z) = x{{z) + ^ mf_,Sf_,{z), 



(1) 



where the principal components, are the eigen- 

functions of the Fisher matrix that describes the depen- 
dence of the polarization spectra on Xe{z) (again, see 
Ref. [l3|), 'm-fi are the amplitudes of the principal com- 
ponents for a particular reionization history, and xl{z) 
is the WMAP fiducial model at which the Fisher matrix 
is computed and from which the principal components 
are obtained. In what follows we use the publicly avail- 
able Sfj_{z) functions and varied the amplitudes for 
ji — 1, 5 for the first five eigenfunctions. Hereafter we 
refer to this method as the MH (Mortonson-Hu) case. 

In the second approach, we use a different parametriza- 
tion, sampling the evolution of the ionization fraction 
as a function of the redshift z at 7 points {z — 
9,12,15,18,21,24,27), and interpolating the value of 
Xe{z) between them with a cubic spline. For 30 < z 
we fix Xe = 2 X 10~^ as the value of Xe expected before 
reionization (and after primordial recombination), while 
Xe = 1 for 3 < z < 6 and Xe = 1.08 for z < 3 in order to 
be in agreement with both Helium ionization and Gunn- 
Peterson test observations. This approach is very similar 
to the one used in Ref. [l^ , and we will refer to it as the 
LWB (Lewis- Weller-Battye) case. 

Exploiting two different and more general reionization 
schemes will provide us with not only a consistency check 
of the results, but also a broader set of reionization his- 
tories. Finally, we also consider here the standard reion- 
ization model, i.e., sudden reionization. 

We have then modified the Boltzmann CAMB code 
[loj incorporating the two generalized reionization sce- 



narios and extracted cosmological parameters from cur- 
rent data using a Monte Carlo Markov Chain (MCMC) 
analysis based on the publicly available MCMC package 
cosmomc [23|. 

We consider here a flat ACDM universe described by 
a set of cosmological parameters 



{tJ6,Wc,e3,n,log[10^"A,]} 



(2) 



where Wb = Q,i,h^ and lUc = ^ch^ are the physical baryon 
and cold dark matter densities relative to the critical den- 
sity, s is the ratio of the sound horizon to the angular 
diameter distance at decoupling. As is the amplitude of 
the primordial spectrum, and n is the scalar spectral in- 
dex. In one case we will also consider the possibility of a 
redshift-dependent dark energy component evolving with 
equation of state w{z) = wq -\- Wzz/{1 + z). 

The extra parameters needed to describe the reioniza- 
tion are the five amplitudes of the eigenfunctions for the 
MH case, or the five amplitudes in the five bins for the 
LWB case, and one single common parameter, the optical 
depth, T, for the sudden reionization case. 

Our basic data set is the seven-year WMAP data 
Q (temperature and polarization) with the routine 
for computing the likelihood supplied by the WMAP 
team. Together with the WMAP data, we also aug- 
ment the WMAP7 data with the CMB datasets from 
BOOMERanG ^i], QUAD d, ACBAR [6], and BICEP 
0] . For all these experiments we marginalize over a pos- 
sible contamination from the Sunyaev-Zeldovich compo- 
nent, rescaling the WMAP template at the correspond- 
ing experimental frequencies. We therefore consider two 
cases: we first analyze the WMAP data alone, referring 
to it as to the "WMAP 7" case, and we then include 
the remaining CMB experiments ("CMB All"). Finally 
we also consider, separately, the galaxy clustering results 
from LRG-7 of Ref. 2l| and the baryonic acoustic oscil- 
lation (BAO) data from the same dataset, see Ref. [22| . 

Our main results are reported in Table I, where we 
list the constraints on n for the different reionization sce- 
narios using different datasets. Notice that using the 
MH parametrization drastically alters the constraints on 
n. For the WMAP7 case alone, the HZ spectrum is not 
only in agreement with the data, but is also close to the 
best-fit value. When the remaining CMB experiments 
are included, the best-fit value of n shifts to lower val- 
ues, but it is still consistent at better than 95% c.l. with 
n = 1. When information from galaxy clustering is in- 
cluded, a value of n = 1 is excluded at about the 95% 
C.L, i.e., in a much less conclusive way than in the sud- 
den case where n = 1 is excluded at more than 99.7% c.l. 
(see e.g., Ref. [isj). 

Adopting the LWB parameterization of reionization 
has very similar effects. The HZ spectrum is only 
marginally excluded at a 95% confidence level if the com- 
plete CMB dataset is considered. Notice that in the 
case of the LWB parametrization, the results are closer 
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FIG. 1: Contours of the 68% and 95% c.l. constraints in the n 
vs T plane for different datasets. The contours regions come 
from a generahzed reionization scenario using (from right to 
left) the MH approach from WMAP-7 data (Yellow), "CMB 
all" (Green), All CMB+BAO (Blue). 



to those obtained in the standard sudden reionization 
scheme. Therefore, the LWB parameterization is less 
general and samples fewer reionization scenarios than the 
MH method. 



Dataset 


Ionization 


n (68% c.l.) 


95% c.l. 


WMAP7 


sudden 


0.965 ± 0.014 


n < 0.993 


CMB All 


sudden 


0.959 ± 0.013 


n < 0.984 


WMAP7 


MH 


0.993 ± 0.023 


n < 1.042 


CMB All 


MH 


0.975 ±0.017 


n < 1.011 


CMB All+LRG-7 


MH 


0.966 ± 0.014 


n < 0.994 


CMB All+BAO 


MH 


0.965 ± 0.014 


n < 0.995 


CMB All+BAO 


MH+u'(2) 


0.985 ± 0.018 


n < 1.025 


WMAP7 


LWB 


0.977 ±0.018 


n < 1.01 


CMB All 


LWB 


0.963 ±0.015 


n < 0.998 



slightly higher than in the standard analysis but consis- 
tent with several physical reionization models. Notice 
also that the HZ spectrum (n = 1) is perfectly com- 
patible with WMAP7. Including more datasets shrinks 
and shifts the contours towards lower values of n, but 
the HZ case is always reasonably consistent at the 95% 
C.L. Including BAO data rules out the HZ spectrum 
at about 2 standard deviations that is considerably less 
stringent than in the standard case. Moreover, consid- 
ering a redshift-dependent equation of state weakens the 
geometrical probes and HZ is again inside the 68% con- 
fidence level for the CMB+BAO case. 

It is important to investigate if a more general reioniza- 
tion scenario with the assumption of n = 1 could provide 
a viable cosmology, i.e., if the value of cosmological pa- 
rameters such as the baryon density or the age of the 
universe are compatible with complementary cosmolog- 
ical information from Big Bang Nucleosythesis and age 
constraints. It has been shown, for example, in Ref. [l3| 
that the assumption of an HZ spectrum in the standard 
analysis with sudden reionization gives a baryon abun- 
dance that is at odds with current bounds from BBN. 
We have therefore performed an analysis using WMAP7 
and assuming n ~ 1 for the MH reionization scenario, 
and the results for the different cosmological parameters 
are summarized in Table 2. 



Parameter 



Constraint (68% c.l.) 



Ho (km s-^Mpc-i) 
Age (Gyr) 



0.0234 ± 0.0004 
0.106 ±0.005 
74.2 ±2.1 
13.6 ±0.1 
0.765 ± 0.022 



TABLE II: Constraints on cosmological parameters from 
WMAP7 assuming a HZ primordial spectrum and a reion- 
ization scenario characterized by the MH parametrization. 



TABLE I: Constraints and upper limits on n from several 
datasets in different recombination scenarios. 



Considering the WMAP7 data alone, the HZ spec- 
trum is at A(— 21n(£)) = 0.72 with respect to the best 
fit model in the case of a MH recombination scenario 
(A(— 21n(£)) = 2.2 when all CMB experiments are con- 
sidered), compared with A(— 21n(£)) = 3.3 when recom- 
bination is sudden (A(— 21n(£)) = 5.8 in case of CMB 
All). 

It is interesting to consider the constraints on the op- 
tical depth T, derived by integrating Xe{z) up to z = 32. 
Figure 1 shows the 68% and 95% c.l. constraints in the n 
vs. T plane arising from different datasets and assum- 
ing the MH reionization parametrization. As we can 
see the optical depth is always in the range 0.08-0.11, 



Notice from the results of Table 2 that the constraints 
on the remaining cosmological parameters are all viable 
when compared with independent cosmological observ- 
ables. Current Big Bang Nucleosynthesis constraints 
from observations of Deuterium abundances provide a 
constraint of Qbh'^ = 0.0213 ± 0.002 at 95% c.l. [H], 
which is compatible with the WMAP7 constraint assum- 
ing HZ. The constraints shown in Table 2 on Hq are also 
well in agreement with the recent HST determination 
h = 0.748 ± 0.036 at 68% c.l. [H. Figure 2 depicts the 
constraints from WMAP7 in the vs. Hq plane in 

a generalized recombination scenario, assuming both a 
varying n and a fixed n — 1 scenarios. Notice that fixing 
n — 1 implies larger values of both ^},h^ and Hq. 

In conclusion, the details of the reionization processes 
in the late universe are not very well known. In the ab- 
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found that Planck will be able to discriminate it from 
a HZ spectrum at more than five standard deviations 
even in the generalized reionization cases describes here. 
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FIG. 2: 68% and 95% c.l. WMAP7 constraints in the Ho vs. 
n^/i^ plane in a generalized reionization scenario described by 
the MH parametrization. The larger (blue) contours regions 
arise when n is allowed to vary. The smaller (red) contours 
correspond to the HZ (n — 1) spectrum. 



sence of a precise, full redshift evolution description of 
the ionization fraction during the reionization period, a 
simple parametrization, with a single parameter z^, has 
become the standard reionization scheme in numerical 
analyses. However, more general reionization scenarios 
are certainly plausible and their impact on the cosmo- 
logical constraints should be carefully explored. For ex- 
ample, processes of dark matter annihilation as in Ref. 
[25| could easily be present. 

In this paper we have investigated the stability of the 
CMB constraints on n in generalized reionization scenar- 
ios. Our study is motivated by the recent claim that 
an HZ spectrum is ruled out by current WMAP data at 
more than 95% c.l. In agreement with previous studies, 
we have found that a more general treatment of reion- 
ization drastically weakens current CMB constraints on 
n. If only the WMAP 7-year data is considered, a HZ 
n = 1 spectrum lies well within the 68% c.l. allowed re- 
gion. If all current CMB datasets are considered in the 
analysis, the HZ spectrum still lies well within the 95% 
c.l. allowed region. If the assumption of constant dark 
energy is relaxed, HZ is again inside the 68% c.l. There- 
fore, current data do not yet rule out an HZ spectrum in 
a conclusive way. It is also worth mentioning that non 
standard processes during recombination (see e.g., Ref. 
[2^ ) or extra relativistic particles (see e.g., Ref. 
could be present and could also put ?i = 1 in agreement 
with current observations. 

Near future data from the Planck satellite mission are 
clearly needed to solve this important question. In this 
respect, performing a parameter analysis with simulated 
data with the experimental configuration described in 
Ref. ^28], and assuming a model with n ~ 0.965, we 
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